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ABSTRACT 

Measurements of the SNe la Hubble diagram which suggest that the universe is ac- 
celerating due to the effect of dark energy may be biased because we are located in 
a 200-300 Mpc underdense 'void' which is expanding 20-30% faster than the average 
rate. With the smaller global Hubble parameter, the WMAP-5 data on cosmic mi- 
crowave background anisotropies can be fitted without requiring dark energy if there 
is some excess power in the spectrum of primordial perturbations on 100 Mpc scales. 
The SDSS data on galaxy clustering can also be fitted if there is a small component 
of hot dark matter in the form of 0.5 eV mass neutrinos. We show however that if the 
primordial fluctuations are gaussian, the expected variance of the Hubble parameter 
and the matter density are far too small to allow such a large local void. Nevertheless 
many such large voids have been identified in the SDSS LRG survey in a search for the 
late-ISW effect due to dark energy. The observed CMB temperature decrements im- 
ply that they are nearly empty, thus these real voids too are in gross conflict with the 
concordance ACDM model. The recently observed high peculiar velocity flow presents 
another challenge for the model. Therefore whether a large local void actually exists 
must be tested through observations and cannot be dismissed a priori. 

Key words: cosmic microwave background, cosmological parameters, cosmology: 
theory, dark matter, large-scale structure of Universe 



1 INTRODUCTION 

The Einstein-de Sitter (E-deS) universe with fi m = 1 is 
the simplest model consistent with the spatial flatness ex- 
pectation of inflationary cosmology. However, Type la su- 
pernovae (SNe la) at redshift z ~ 0.5 app ear ~ 25% 
fainte r than expected i n an E-deS universe IjRiess et al.l 
1 19981 ; iPerlmutter et all 1 19991 ). Together with measure- 
men ts of galaxy clustering in the Two-degree Field sur- 
vey ( Efstathiou et al. 2002) and of cosmic microwave back- 
ground (CMB) anisotropies by the Wilkinson M icrowave 
Anisotropy Probe (WMAP) (|Spergel et al.ll2003l ). this has 
established an accelerating universe with a dominant cos- 
mological constant term (or other form of 'dark energy') 
which presumably reflects the present microphysical vac- 
uum state. This 'concordance' ACDM cosmology (with 
£7a — 0.7, Q m ~ 0.3, h ~ 0.7) has passed a number of 
cosmological tests, inc luding baryonic acoustic oscillations 
IjEisenstein et al.l|2005l ) and measurements of mass fluctu- 
ation s from clusters and weak lensing (e.g. IContaldi et al 
2003). Further observations of both SNe la i Riess et al 



12004 lAstier et all |200d IWood-Vasey et ail 120071 ) and the 
WMAP 3-year results i|Spergel et al.l 120071 ) have continued 
to firm up the model. However there is no physical basis for 
this model, in particular there are two fundamental prob- 
lems with the notion that the universe is dominated by vac- 
uum energy. The first is the notorious fine-tuning problem 
of vacuum fluctuations in quantum field theory — the en- 
ergy scale of the cosmological energy density is ~ 10~ 12 
GeV, many orders of magnitude below the energy scale of 
~ 10 2 GeV of the Standard Model of particle p hysics, not 
to m ention the Planck scale of ~ 10 19 GeV (see IWeinberd 
Il999f ). The second is the equally acute coincidence problem: 
since pa/Piu evolves as the cube of the cosmic scale factor 
a, there is no reason to expect it to be of 0(1) today, yet 
this is apparently the case. In fact what is actually inferred 
from observations is not an energy density, just a value of 
O(Hq) for the otherwise unconstrained A term in the Fricd- 
mann equation. It has been suggested that this may simply 
be an artifact of interpreting cosmological data in the (over- 
simplified) framework of a perfectly homogeneous universe 



in which H ~ 10~ 42 GeV - 
the problem (|Sarkarll200sT ). 



(10 cm) is the only scale in 
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ergy if the assumption of a scale-invariant primordial power 
spectrum is relaxed. This assumption is worth examining 
given our present ignor ance of the physics behind inflation. 
We have demonstrated IjHunt fc Sa rkar 2007) that the tem- 
perature angular power spectrum of an E-deS universe with 
h ~ 0.44 matches the WMAP data well if the primordial 
power is enhanced by ~ 30% in the region of the second 
and third acoustic peaks (corresponding to spatial scales 
of k ~ O.Of - 0.1 h Mpc" 1 ). This alternative model with 
no dark energy actually has a slightly better \ 2 f°r the fit 
to WMAP-3 data than the 'concordance power-law ACDM 
model' and, inspite of having more parameters, has an equal 
value of the Akaike information criterion used in model se- 
lectio n. Other E-deS models with a broken power-law spec- 
trum (JBlanchard et al.ll2003h have also been shown to fit the 
WMAP data. Moreover, an E-deS universe can fit measure- 
ments of the galaxy power spectrum if it includes a ~ 10% 
component of hot dar k matter in the form of massive neutri- 
nos o f mass ~ 0.5 eV l|Hunt fc Sarkarll2007l ; lBlanchard et al.l 
120031 ) . Clearly the main evidence for dark energy comes from 
the SNe la Hubble diagram. 

A mechanism that sets A = is arguably more plausible 
than one which leads to the tiny energy density pA — 10 -47 
GeV 4 associated with the concordance cosmology|j If A is in- 
deed zero then perhaps some effect fools us into wrongly de- 
ducing the existence of dark energy by mimicking a nonzero 
cosmological constant. It is natural to connect this effect 
with inhomogenities since cosmic acceleration and large 
scale nonlinear structure formation appear to have com- 
menced simultaneously. This approach offers the possibil- 
ity of solving the cosmological constant problems within the 
framework of general relativity and keeps the introduction 
of new physics to a minimuml 2 ] Several different ways in 
which inhomogenities could potentially mimic dark energy 
have be e n con si dered in the li t erature — for reviews see 
ICelerierl l|2007h ; iBuchertl [[20081 ); lEnqvistl IJ2008D . In an in- 
homogeneous universe averaged quantities satisfy modified 
Friedmann equations which contain extra terms correspond- 
ing to 'backreaction' since the operations of spatial averag- 
ing and time evolution do not commute (Buchcrt 2000). The 
backreaction terms depend upon the variance of the local ex- 
pansion rate and hence increase as inhomogenities develop. 
Whether backreaction can indeed account for the apparent 
cosmological acceleration is hotly debated and remains an 
open question at presen t (Wcttcrich 2003; Ishibashi & Wald 



2006; Vandcrvcld ct al. 2007; Wiltshire 2007; Kh osravi et al. 
20071: iLeith et al.ll2008l; iBehrend etaljbooa : iRasanenlbo'ol 
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Another possibility is that inhomogeneities affect light 
propagation on large scales and cause the luminosity 
distance-redshift relation to resemble that expected for 
an accelerating universe. This has been investigated for a 
'Swiss-cheese' universe in which voids modelled by patches 



1 'Quintessence' models, which attempt to address the coinci- 
dence problem, also assume that every other contribution to the 
vacuum energy cancels apart from that of the quintessence field. 

2 In models that seek to explain the observations through modifi- 
cations of gravity, the relevant scale of Hq has to be introduced 
by hand, just as in quintessence models the quintessence field has 
to be given a mass of order Ho — these are technically unnatural 
choices since this is an infrared scale for any microphysical theory. 



of Lemaitre-Tolman-Bondi (LTB) space-time are distributed 
throughout a homogenous background. However, the results 
seem to be model dependent: some authors find the change 
in light pro pagation to be negligib l e because of cancella- 
tion effects jBiswas fc Notaril [20081 : IBrouzakis et alj |20Q; ' 
iBrouzakis fc Tetradisl l200Sf ). whereas iMarra et al.l (| 20o" 
claim it can par tly mimic dark en e rgy if the voids have 
radius 250 Mpc (|Marra et all 120081 '). iMattssonl |2007f ) has 
noted that observers may preferentially choose sky regions 
with underdense foregrounds when studing distant objects 
such as SNe la, so the expansion rate along the line-of-sight 
is then greater than average; such a selection effect he argues 
can allow an inhomogeneous universe to fit the observations 
without dark energy. 

In this paper we are mainly interested in a 'local void' 
(sonretimes referred to as "Hubble bubble") as an explana- 
tion for dark energy; to prevent an excessive CMB dipole 
moment due to our peculiar velocity we must be located 
near the centre of the void. An underdense void expands 
faster than its surroundings, thus younger supernovae in- 
side the void would be observed to be receding more rapidly 
than older supernovae outside the void. Under the assump- 
tion of homogeneity this would lead to the mistaken conclu- 
sion that the expansion rate of the Universe is accelerating, 
although both the void and the global universe are actu- 
ally decelerating. Henceforth we use the 'Hubble contrast' 
8h = (H[ n — H out ) /H out to characterise the void expansion 
rate, where Hi n and -ff ou t are the Hubble parameters inside 
and outside the void respectively. (Other authors have used 
the 'jump' J = Hi n /H ou t — 1 + 8h to characterise the 
void.) The reduced Hubble parameter h is defined as usual 
by -ff ou t = 100ft km s _1 Mpc" 1 throughout. 

The local void scenario has been investig ated by 
several authors us i ng a vari ety of methods llCeleriei 




2001a lblfcT; [iguchi et al] 120021: 



lAlnes et alj 



Tomita 
Mansouri et all 120011; 



Vanderveld et all I2006J : ICarfinkld [20061: IChung fc Romano 



Alexander ct al. 



12006 

2007; 



Alnes fc Amarzg uioui 

Biswas et all 120071: 

Caldwell fc Stebbinsll2008l: IClarkson et a l. 2008; Uzan et al. 

Garcia-Bellido fc Haugboelld l2008al : IClifton et ~ 



Garcia-Bellido fc Haugboelld l2008bh . In a series of 
papers, Tomita modelled the void as a open Friedmann- 
Robertson- Walker (FRW) region joined by a singular mass 
shell to a FRW background and found that a void with 
radius 200 Mpc and Sh = 0.2 5 fits the supe r nova Hubble 
diagram without dark energy (lTomitall200ld ). lAlnes et al.l 
(2006) showed that a LTB region which reduces to a 
E-deS cosmology with h — 0.51 at a radius of 1.4 Gpc 
with 8h = 0.27 can match both the supernova data 
and the location of th e first acoustic peak in the CMB. 
lAlexander et al] ([2007M attempted to find the smallest 
possible void consistent with the current supernova results 
— their LTB-based 'minimal void' model has a radius of 
350 Mpc and J ~ 1.2 i.e. Sh — 0.2; a void o f similiar size 
but w ith Sh — 0.3 had been discussed earlier l|Biswas et al] 
[2007J). Unfortunately, since this model is equivalent to an 
E-deS universe with h = 0.44 outside the void where the 
Sloan Digital Sky Survey (SDSS) luminous red galaxies 
lie, as it stands it is unable to fit the measurements of 
the baryonic acoustic oscillation (BAO) peak at z ~ 0.35 
l|Blanchard et al.ll2006T ). LTB models of much larger voids 
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were considered by iGarcia-Bellido fe Haugboelld (|2008al ) 
(with radii of 2.3 Gpc and 2.5 Gpc and Hubble contrasts of 
0.18 and 0.30 respectively) and it was demonstrated they 
can fit the supernova da ta, BAO data and the location 
of the first CMB peak. IClifton et all l|200Sl ) found the 
best fit to the SNe la data for a void of radius 1.3 ± 0.2 
Gpc and an underden s ity o f about 70% at the centre 
and Bolci ko fe Wvithe I (J2008J) confirmed that such a void 
provides an excellen t fit to th e latest 'Un i on da taset' 
JKowalski et~aTll2008l ). Moreover linoue fe Silkl (|2006l ) have 
shown that the unexpected alignment of the low multipoles 
in the CMB anisotropy can be attributed to the existence 
of a local void of radius 300 ft" 1 Mpc. These authors also 
suggested that the anomalous 'cold spot' in the WMAP 
southern sky is due to a similar void at z ~ 1 and some 
evide nce for this has emerged subsequently iJRudnick et al.l 
120071 ). Recently, a large number of voids of varying sizes 
have been identified in the SDSS Luminous Red Galaxy 
(LRG) catalog in a search for the late int egrated Sachs- 
Wolfe (ISW) effect due to dark energy (JGranett et al.l 
l2008bh . 

How likely is the existence of such huge voids ac- 
cording to standard theories of structure formation? 
Statistical measures of the void distribution such as 
the void probability function and underdense probabil- 
ity function have been estimated from th e 2dfGRS, SPSS 
and DEEP2 galaxy redshift surveys IIHovle fc Vogelev 



2004 ; ICroton et all |2004 IPatiri et all 120051; IConroy et al 



2003: ITikhonovl 120061: iTinker et all 120071 : iTikhonovl l2007t 
von B enda-Bec kmann fc Mueller! |2007). Void probability 
statistics have also be en examined theoretically us- 
ing analytical methods (ISheth fc van de Weveaertl 120041 ; 
iFurlanetto fc Piranl liooH; IShandarin et all 120061) and N- 



, , 



body simulations (iLittle & Weinberg 1994; Sch midt et all 
20001; lArbabi-Bidgoli fc Muellerl 120021; iBenson et all booa 
Padilla et all 120051 ). However such studies have been re- 



stricted to voids with radii of 10-30 Mpc. The scales of the 
large voids we are considering lie in the linear regime where 
the variance of the Hubble contrast is directly related to 
the matter powe r spectrum ? m (fc). It has been noted (us- 
ing results from iTurner et all (119921 )') that above 100 Mpc 
linear theory predictions agree well with N-body simulation 
results, although on smaller scales the Hubble contrast is un- 
derestimated due to non-linear effects (|Shi et al.ll 19961 ) . Ap- 
plying ljrwar_theory andusing the measured CMB dipole ve- 
locity, IWang et all (|l998f ) obtained the model-independent 
result {Sh}j( < 10.5 h^ 1 Mpc/R in a sphere of radius R. 
(This ought to be an acceptable procedure up to scales of or- 
der 800 /i _1 Mpc — on larger scales, relativistic corrections 
become increasingly important.) In this paper we update 
these results by determining the probability distribution of 
8h and the density contrast on various sc ales using con- 
straints on P m (fc) from WMAP 5-year data (Komatsu et al 



2008) and the SDSS galaxy power spectrum JTegmark et al 
2003). We find that even the 'minimal local void' is ex- 



tremely unlikely if the primordial density perturbation is 
indeed gaussian as is usually assumed and the other LTB 
model voids even less so. However by the same token, the 
ISW effect due to the voids seen in the SDSS LRG survey 
JGranett et al.ll2008bl ) appears to be too strong. Moreover, 
observed large-scale peculiar velocities appear to be much 
higher than expected JKashlinskv et alll2008l ; IWatkins et all 



2008). It would appear that the standard model of structure 
formation itself needs reexamination hence the existence of 
a large local void cannot be dismissed on these grounds. 



2 MODELS 

We study variations of the Hubble parameter in the context 
of two different cosmological models, both of which fit the 
WMAP and SDSS data but have different amounts of power 
on spatial scales of 0(100) Mpc. The intention is to examine 
whether previou s conclusions conc erning the magnitude of 
such variations l|Wang et all 1 1998) can be circumvented in 
an unorthodox model. 

Our first model is the standard ACDM concordance 
model with a power-law primordial power spectrum. The 
spectral index and amplitude Viz of the comoving curva- 
ture perturbation spectrum are evaluated at a pivot point of 
fc = 0.05 Mpc" 1 . The second model is dubbed the 'CHDM 
bump model' since it has both cold and hot dark matter 
and a ' bump' in the primord ial spectrum. It was developed 
by us l|Hunt fc Sarkarl 120071 ) based upon the supergravity 
multiple inflation scenario in which 'flat direction' fields un- 
dergo gauge symmetry-breaking phase transition s during in- 
flation triggered by the f all in temperature (JAdams et all 
ll997l ; lHunt fc Sarkarll2004l ). Each flat direction ip has a grav- 
itational strength coupling to the inflaton <f>, giving a con- 
tribution to the potential of the form V C h\(f> 2 ip 2 - The flat 
directions are lifted by supergravity corrections and non- 
renormalisable superpotential terms. Thus when a phase 
transition occurs the flat direction evolves rapidly from the 
origin where it was trapped by thermal effects to the global 
minimum of the potential. Each phase transition changes the 
effective inflaton mass from m| to m| — A (ip) . Since the pri- 
mordial power spectrum is very sensitive to the inflaton mass 
this can introduce features into the spectrum. We showed 
that two flat directions ip\ and ip2 which cause successive 
phase transitions about 2 e-folds apart and create a small 
bump in the power spectrum centred on fc ~ 0.03 ftMpc -1 , 
allow an E-deS mode l with h = 0.44 to fit the WMAP data 
l|Hunt fc Sarkarll2007l ). The effective scalar potential is: 
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Here ti and £2 are the times at which the first and second 
phase transitions begin, Ai and A2 are the couplings between 
<f> and the flat directions, 71 and 72 are the co-efficients of 
the non-renormalisable terms of order m and 712, and Vo 
is a constant which dominates the potential. In the slow- 
roll approximation the height of the bump is Viz , and the 
amplitude of the primordial perturbation spectrum to the 
left and right of the bump is Viz and Viz respectively, 
where 
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are the fractional changes in the inflaton mass- squared due 
to the phase transitions. The bump lies approximately be- 
tween the wavenumbers k\ and fe where fe = feie i - t2 ~ tl > . 
In this paper we set 71 and 72 equal to unity, m 2 — 0.005-ff 2 , 
0o = O.OIA/p, ni = [4 = 3# 2 an d A! = A 2 = # 2 /M g 
throughout as in our earlier work IjHunt fc Sarkarl 12007). 
In fitting to the WMAP-5 data we also consider con- 
tinuous (non-integral) values of n\ and n 2 to determine 
whether a different shape of the 'bump' gives a better fit, 
keeping in mind that its physical origin may be differ- 
ent fr o m multiple inflation llChung et al . 1999; L esgoi 
19991: iKaloper fc Kaphnghatl [20031: lEasther et all 
" 120051 : iGond 120051; lAshoorian fc Krausel 
2008). 



Wan g et al, 
Bean et al.l 



A pure cold dark matter (CDM) model exhibits ex- 
cessive galaxy clustering on small scales. Therefore it is 
necessary to include a hot dark matter (HDM) compo- 
nent which suppresses structure formation below the free- 
streaming scale. We obtain a good match to the shape of 
the SDSS galaxy power spectrum with 3 neutrino species 
of mass ~ 0.5 eV. Hence t he CHDM bump mo del has 
n b ~ 0.1, n„ ~ 0.1, Q c ~ 0.8 (|Hunt fc Sarkarll2007l y 



3 THE DATA SETS 



We fit to the WMAP 5-year iJNolta et al.ll2008l ) temperature- 
temperature (TT), temperature-electric polarisation (TE), 
and electric-electric polarisation (EE) spectra. Compared to 
the WMAP-3 results, the WMAP-5 measurement of the TT 
spectrum is ~ 2.5% higher in the region of the acoustic 
peaks due to the revised beam transfer functions, and the 
third acoustic peak is determined more accurately. Polari- 
sation measurements are improved by the use of data from 
an additional waveband. 

We also fit the linear matter power spectrum T } m (k) to 
the measurement o f the real space galax y power spectrum 
V g (k) in the SDSS (JTegmark et al.ll2003h . 



4 METHOD 

The Hubble contrast Sh smoothed over a sphere of radius 
i? is (|Shi et alJll996h 



s 1 \R I j3 v(y) 
fo(x) = / d y 



Hout |y - x| 



■ W R (y 



(7) 



where v is the peculiar velocity field and Wr is the 'top hat' 
window function, 



W r (k) 



3/(4nR 3 ), 
0, 



|x| < R, 

Ixl > R. 



(8) 



Using linear perturbation theory (|PeeblesNl993l ) it can be 
shown that the var iance of 8h i s related to the matter power 
spectrum as (jWang et al.lll998T ) 



(*«>«= 2^ 



dkk 2 P m (k)W H (kR). 



Here the window function Wh is 
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and the dimensionless linear gr owth rate / for a ACDM uni- 
verse can be approximated by IjLahav et al.lll99ll : iHamiltonl 
200lf| 



/(n m ,n A ) ^fim /7 + 
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Similarly, the variance of the density contrast S 
(pin — Pout) /pout in a sphere of radius R is 



( S ' 2 )r = ^ 



dk k 2 V m (k) W 2 (kR) . 







where the window function 
3 



W (kR) = 



(sin kR — kR cos kR) , 



k s R s 
is the Fourier transform of Wr. 

The variance of the peculiar velocity is given by 

/ Hout 



(v 2 )r = 



2tt 2 



dkV m (k)W 2 (kR). 



(12) 



(13) 



(14) 



Finally we also consider £l ln — 8irGpi n /H 2 n , the ratio 
of the matter density to the critical d ensity as measure d 
locally by an observer inside the void (JWang et al.lll998l ). 
The variance of the perturbation fo = (fiin — ^m) film is 
then 



(15) 
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(16) 



We use the Monte Carlo Markov Chain (MCMC) ap- 
proach to cosmological parameter estimation, which is a 
method for drawing samples from the posterior distribu- 
tion P(rujdata) of the cosmological parameters vj, given 
the data. For a disc ussion of the MCMC likelihood analysis 
see Appendix B of iHunt fc Sarkarl (J2007I ). Given n samples 
vj^' the best estimate for the distribution is 



n 

P(tu|data) ~ -^5° (z 



Ji)\ 



(17) 



3 IHamiltonl J2001T ) emphasized that the power-law exponent is 
4/7 for a hig h density universe, s o we have corrected the previous 
formula from lLahav et al.l ill 99 if ) accordingly. 
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where 8 D is the Dirac delta function. The CMB angular 
power spectrum and the matter power spectru m (corrected 
for n on- linear evolution using the 'Halofit' ([Smith et al,l 
12002 ) procedure) of each model are calculated using a mod- 
i fied version of t he CAMqj cosmological B oltzmann code 
Lewis et allfeOOOl ) following the approach of lHunt fc Sarkarl 



2007). While the temperature of the CMB monopole would 
be affected if we are located near the centre of a spher- 
ically symmetric void, secondary anisotropies due to the 
void are expected to decay rapidly for higher multipoles 
(| Alexander et al.ll2007l ). We therefore neglect the possible ef- 
fects of the void on the angular power spectra since these are 
important only at low multipoles where the cosmic variance 
is large. The integrals for the variances in eqs.([9l [T2lll5[) were 
evaluated numerically. Care was taken to ensure the precise 
values of the integration limits did not affect the results. We 
compute / num erically using the GROwA software package 
llHamiltorfeOOlF ). We use a version of the COSMOMCO package 
i|Lewis fc Bridlell2002l ) modified to include (<5ff)^ /2 , (S 2 ) 1 ^ , 

(v 2 ) and \5q) for 8 values of R as additional derived 
parameters determined from the base Monte Carlo param- 
eters. It follows from eq.{T7| that I-D marginalised distri- 
butions of these quantities for each R value are obtained by 
plotting histograms of the samples. The probability distri- 
bution P (Sh\ data) r of Sh on the scale R given the data 
can be written as 



(18) 



P (S H \ data,) r ■= / P (S H \ vb) r P {vj\ data) dt 
Using eq. (|17[l this is approximated by 



n 

P (5 H \data) R = -J2 P ( 5 h\ 



m )R, 



where 



P(Sh\ 



■G7)R = 



y/^mk 



exp 



H*h) r 



(19) 



(20) 



We calculate the probability distribution P{8\ data) r in 
same way. 

Flat priors are used on the parameters listed in Table 
[T] Here 9 is the ratio of the sound horizon to the angular 
diameter distance to last scattering (multiplied by 100) , r is 
the optical depth (due to reionisation) to the last scattering 
surface, and /„ = fi„/fid is the fraction of dark matter in 
the form of neutrinos, where the total dark matter density is 
fid = fi c + fii/ • We assume the chains have converged when 
the Gelman-Rubin 'R' statistic falls below 1.02. We evaluate 
the sum in eq. (|19[) when post-processing the chains. 



5 RESULTS 

The mean values of the marginalised cosmological param- 
eters together with their 68% confid ence limits are listed 
in Tabled As in our previous work (Hunt & Sarkar 2007) 
we also list the value of the Akaike information criterion 
(AIC) relative to the ACDM power-law model. Recall that 



4 http://camb.info 

5 http://cosmologist.info/cosmomc/ 



the AIC is defined as AIC = -2 In £ max + 2N (|Akaikelll974r ) 
where £ max is the maximum likelihood and N the num- 
ber of parameters. It is a commonly used guide for judging 
whether additional parameters are warranted given the in- 
creased model complexity, and quantifies the compromise 
between improving the fit and adding extra parameters. 

The CHDM 'bump' model with m = 12 and n 2 = 13 
has a x equal to the ACDM power-law model. Allowing m 
and n2 to vary freely further improves the fit to the data 
with the consequence that the CHDM model with m and 
712 continuous is favoured over the ACDM model according 
to the AIC. The primordial power spectrum of the models 
is shown in FigH together with the fit to the WMAP TT 
and TE spectra and the SDSS galaxy power spectrum. 

The uncertainties of the derived parameters are smaller 
compared to those derived from the WMAP 3-year results, 
as would be expected for higher quality data. For example, 
the optical depth due to reionisation for the CHDM model 
with continuous m and 712 has gone from r = 0.0751q'qi2 



to t = 0.0771 



+0.0073 



due to the more accurate polarisation 



measurements. The shape of the 'bump' in the primordial 
power spectrum for the CHDM model with continuous n\ 
and ri2 is slightly changed by the new data. Although the 
quantity In (l0 10 P TC (O) ) is almost unaltered, In (l0 10 P TC (1) ) 
has increased slightly from a 3-year value of 3.429lJJ'Q4g to 
a 5-year value of 3.462 +q ^ because of the increased ampli- 
tude of the TT spectrum for multipoles I > 200. Due to the 
increased height of the third acoustic peak, In (10 10 Ptj' 2 )) 



has increased from 3.091 



+0.071 



to 3.183 



+0.043 
-0.041 



and 10 4 fc 2 



fallen from 585±|| Mpc" 1 to 500±^ Mpc -1 . The increased 
amplitude of the primordial power spectrum on small scales 
has raised a$ from a value of 0.662± '; o l g to 0-700±g:g|| . 
The mean values of the variances (J>h) „> (5/ „, (v/ „ 



n). 



together with their la limits, are plotted in 



and (5; 

FigEl The different variances in the two models can be un- 
derstood with reference to the matter power spectrum. From 
the relativistic Poisson equation, a given density perturba- 
tion leads to a larger curvature perturbation in a higher 
density universe. Since the amplitude of the primordial cur- 
vature perturbation is similar in both models (as can be 
seen from FigUJ the density contrast during the early mat- 
ter dominated era is greater in the ACDM universe than in 
the higher density CHDM universe. Although the growth 
of density perturbations at late times is suppressed in a low 
density universe, this means that the matter power spectrum 
of the ACDM universe is larger on all scales than that of the 
CHDM universe, when measured in units of h~ 3 Mpc 3 . (This 
is not evident in Fig[T] where the galaxy power spectrum is 
shown — the galaxies are more biased in the CHDM universe 
than in ACDM so the matter power spectrum is lower.) 

This also explains why, as seen in Fig(2] (5 2 ) is uni- 
formly greater for the ACDM model. The linear growth fac- 
tor / is smaller for the ACDM universe, and the peak in 
the matter power spectrum occurs at a larger scale. Thus 
the quantity / 2- P m (fc) which appears in eq.® is greater 
for the ACDM universe for wavenumbers below fc cr oss — 
0.01 /iMpc -1 but is greater for the CHDM universe for 
wavenumbers above k CIOSS . The window function Wh (|10l) 
makes (5%) sensitive to the value of f 2 Vm (fc) for the 
wavenumber fc ~ n/R. Consequently the \S%) curves for 
the two models cross at the scale 7r/fc cross ~ 300 ft -1 Mpc 
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Table 1. The priors adopted on the base Monte Carlo parameters of the various models, as well as on the derived 
parameters: the Hubble constant and the age of the Universe. 



Parameter 



Model 



ACDM power-law 



CHDM bump with 
ni = 12, n,2 = 13 



CHDM bump with 
n±, ri2 continuous 





Lower limit 


Upper limit 


Lower limit 


Upper limit 


Lower limit 


Upper limit 


n b h 2 


0.005 


0.1 


0.005 


0.1 


0.005 


0.1 


n c h 2 


0.01 


0.99 










8 


0.5 


10.0 


0.5 


10.0 


0.5 


10.0 


T 


0.01 


0.8 


0.01 


0.8 


0.01 


0.8 


fu 






0.01 


0.3 


0.01 


0.3 


n 3 


0.5 


1.5 










I0 4 k 1 /Mpc- 1 






0.01 


600 


0.01 


600 


10 4 fc 2 /Mpc- 1 






0.01 


800 


0.01 


1100 


In (10 10 Vk) 


2.7 


1.0 










\n(l0 10 Vn (0) ) 






2.0 


6.0 


2.0 


6.0 


ln(l0 10 Vn {1) ) 










2.0 


6.0 


ln(l0 10 7V 2 >) 










2.0 


6.0 


h 


0.4 


1.0 


0.1 


1.0 


0.1 


1.0 


Age/Gyr 


10.0 


20.0 


10.0 


20.0 


10.0 


20.0 



as seen in Fig[2] The two (v ) curves cross at a smaller 
scale of about 100 ft -1 Mpc. This is because the integral 
(|14[) for the variance in the peculiar velocity (v 2 ) is more 
strongly weighted towards small wavenumbers than the cor- 
responding expression eq.© for the variance in the Hubble 
contrast (S H ) , which has an additional factor of k 2 . Fi- 
nally for \Sn)„ the situation is intermediate between that 
for the variance in the density contrast and the variance in 
the Hubble contrast, since only some of the terms in Wn 
contain factors of /. 

The scale dependence of \5 H / R is the reason that the 
P (Sh I data) distribution is broader for the ACDM power- 
law and the CHDM 'bump' models on scales above and be- 
low 300 h~ x Mpc, respectively, as shown in Fig(3] Similarly 
the P (5\ data)_R distribution is broader for the ACDM model 
on all scales, as seen in Fig (4] 

To illustrate our findings we calculate the probability of 
a fluctuation in the Hubble contrast greater than or equal 
to a given value 8 H in a sphere of radius R, given by: 



Probability (Sh ^ S H 



model on small scales because the P (Sh | data)/? distribu- 
tion is broader for the CHDM model on these scales. Con- 
versely since the distribution is broader on large scales for 
the ACDM model, the probability is greater there for this 
model. 

Similarly we calculate the probability of a fluctuation 
in the density contrast less than or equal to a given value 
—5° in a sphere of radius R, which is given by: 



Probability (S < S°) R = / P (8\ data) H d<5. 

J — oc 



(22) 



This probability is greater for the ACDM model on all scales 
as seen in Fig(S] due to the broader P (S\ data)_R distribution. 
Moreover, we can determine the probability of one or 
more voids with comoving volume Vi occurring within some 
larger comoving volume V2. If the ratio V2/V1 is N to the 
nearest integer and p is the probability of a void with 
volume Vi, then the probability of n voids within V2 is 
( n )p n (1 — p) " where (j is the binomial coefficient. The 
expected number of voids within V2 is Np. 



P (S H \ data) r dSn- 



(21) 



Since P (5h\ data)_R is symmetric this is also equivalent to 
the probability of a fluctuation being less than or equal to 
—5 H . As seen in FigO the probability of a large excursion 
in Sh is largest on small scales, in accordance with physical 
intuition. Note that the probability on all scales tends to 
a value of 1/2 for small S H because the fluctuation has an 
equal probability of being positive or negative. The proba- 
bility is greater for the CHDM model than for the ACDM 



6 DISCUSSION 

A void with Sh — 0.2 — 0.3 and a radius exceed- 
ing 100 /i _1 Mpc is r equired to fit t he supernova data 
without dark energy |Tomital l200ld : iBiswas et ail 120071 ; 
I Alexander et al.ll2007l ). The probability that we are situated 
in such a void is less than 10 -12 as can be seen from Fig[5] 
The probability is exponentially smaller for the larger voids 
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Table 2. The marginalised cosmological parameters for the various models (with lcr limits). 
The 12 parameters in the upper section of the Table are varied by CosmoMC, while those in 
the lower section are derived quantities. The \ 2 °f the fit is given, as is the Akaike information 
criterion relative to the power-law ACDM model. 



Parameter 



ACDM power-law 



Model 
CHDM bump with CHDM bump with 
n± = 12, ?i2 = 13 Hi, n<z continuous 



n h h 2 






n n99^4+oooo 60 

U.UZiO*_ .00061 


01674+ - 00041 

U.U1D/1_ 000 4 7 


01762+ ' 00095 


Q c h 2 






1U4+ - 0046 
u.ii44_ 0046 






e 






1 0397+ - 0029 
1 ' uot " -0.0031 


i rmi+O' 0039 

1.UO±1_q 003 g 


i rww+O- 0048 
i.uoaz_ 0047 


T 






n n«49+ - 0077 
0.0842_ 00082 


0721+ 00069 
u.u(zi_ 0Q75 


0771+ - 0073 
U.U( li_ .0083 


u 








H4+0-015 
U - ±J ^-0.012 


0.085_ ; 022 


n B 






Q61+ ' 014 

U.MOl_ Q14 






10 4 fci/Mpc~ 


1 




8i.7t|;« 


vr±\\ 


10 4 fc 2 /Mpc~ 


1 




442±g 


500t 2 4 


ln(lO : 


i0 Tv) 




o n7 o+0.037 
°' UIO -0.037 






ln(lO : 


l0 Vn m ) 




, 2g4 +0.031 
°""' -0.031 


o 974 +0.048 
°- z,y: -0.048 


In (l0 : 


w Vk (i 


J ) 






3 462+0.036 
o.tvz_ Q 036 


ln(lO : 


10^(3 









o ion+0.043 
o.ioo_ 041 


Q c h 2 








1450+ 00079 
u.iiou_ 0077 


U.10D_ 013 


Q d h 2 








1634+ ' 0042 

U.1DJ4_ 0Q45 


1702+ - 0073 


h 






0.695±°;°£ 


4244+0.0052 
u -* z4t4 -0 .0055 


fl 4QQQ+0.0093 
U -* MO -0.0094 


Age/Gyr 




13.78+^ 


15 S6+ - 20 
lO-^D-0.19 


15.05_ 32 


i I iYi 






n 9S/I+0.025 






n A 






716+ ' 025 
u ' ,1D -0.025 






OS 






o.8i7±g;g? 


617+ - 059 
u - Di '-0.055 


700+ 0098 
u.(uu_ ogg 


^rcion 






n-o±i:l 


13 0+ 20 
io.u_ 2 


134+21 


Am 2 








07495+ - 00046 
U.U/*»O_ 00046 


n n«Q+ 002 o 


Am 2 , 








n 1 51 QQ+0. 00084 
U.1J1DO_ ,00084 


136+ 0015 

U.1JD_ 016 


H(t 2 - 


-*l) 






i-68i° :g 


■t 7 o + 0.14 
J ~'°-0.18 


x 2 






1339.9 


1339.9 


1330.2 


Aaic 









6.0 


-3.7 



of Gpc size that have also been considered IIAlnes e t al. 2006; 
iGarcia-Bellido fc Haugboelldl2008al ; IClifton et al.ll2008l )l 6 l 

However before we dismiss the possibility of a local void 
on these grounds we should also evaluate the probability of 



6 There is a further constraint on Gpc scale voids from the 
observ ed absence of a 'y-distor tion' in the spectrum of the 
CMB dCaldwell fc Stebbinsll2008h and from the 'kinetic Sunyaev- 
Zeldovich' effect obs erved for X-ra y emitting galaxy clusters 
IIGarcia-Bellido fc Haugboel lc 2008b). However this has no im- 
pact on smaller voids. 



voids which have actually been claimed to exist elsewhere 
in the universe. For example it has been argued that a void 
with radius 200 — 300 h~ Mpc and an density contrast of 
8 = —0.3 at z ~ 1 c an account for the W MAP 'cold spot' in 
a ACDM universe i|lnoue fc S ilk 2006). Even if we conser- 
vatively take the radius to be 150 /i _1 Mpc (and the same 
underdensity) , the probability that one or more such voids 
lie within the volume out to z = 1 is only I.O5I093 x 10 -10 - 

It has been argued t hat the WMAP cold spot may 
not be a localized feature ([Naselskv et al.l 120071 ) and there 
may be no matching void in the NVSS radio source cat- 



8 Paul Hunt and Subir Sarkar 



3.2 




oonor 



^ 4000- 



21)00- 



CHDM bump 
ACDM power-law 




io-" io- 

/.-(/iMpcT 1 ) 



10 100 

Multipole moment (/) 



1000 



1.0 



b 0.5 
- 1 



-0.5 




- CHDM bump 
— ACDM power-law 



H. 10 J 





10 100 

Multipole moment (/) 



1000 10 



- CHDM bump 
— ACDM power-law 

J^/iMpc -1 ) 




10 



Figure 1. The top left panel shows the primordial perturbation spectrum for the CHDM bump model (with n\ = 12 
and 712 = 13) and for the ACDM power-law model with n B ~ 0.96. The top right and bottom left panels show the 
best-fits for both models to the WMAP-5 TT and TE spectra, while the bottom right panel shows the best-fits to 
the SDSS galaxy power spectrum. 



alogue IjSmith fc Hutererll2008r ). however an equally strik- 
ing anomaly arises if we consider the large number of voids 
which have been identified i n the SDSS LRG su rvey in a 
search for the late ISW effect (|Granett e t al. 2008a b). These 
are of angular radius ~ 4° corresponding to a (comoving) 
radius of ~ 50 h^ 1 Mpc and are tabulated as having la, 
2a or 3a underdensities. These numbers relate to the de- 
tection significance (the likelihood of detecting the void by 
chance out of a Poisson distribution) rather than the likeli- 
hood of finding such underdensities in a gaussian field which 
we have computed in this paper (B. Granett, private com- 
munication). Moreover the observed LRGs are biased with 
regard to the dark matter hence the underdensities in dark 
matter are likely to be smaller t han th e quoted values. 

However if iGranett et al,l i|2008al lbl) have indeed de- 
tected the late ISW effect as they assert, we can simply 
circumvent these uncertainties by requiring that the voids 
be large enough and/or underdense enough to yield the ob- 
served CMB temperature decrements. To calculate the late 
ISW effect we consider the propagation of CMB photons to 
us from the last scattering surface through an intervening 
void. The photon temperature change caused by the void is 



AT 



d$ , 
— da, 

da 



(23) 



where af ar is the scale factor when the photon crossed the 
far side of the void and a ncar is the scale factor when the 



photon crossed the near side of the void. The gravitational 
potential of a void with proper radius r is 

4ttG 



$ 



-r p h 8 (a) . 



(24) 



Here the background density is given by pb = 
3//of2 m /87rGa 3 and the density perturbation is given by 
8 (a) — D (a) 8 (ao) where D is the linear growth factor. 
Hence 

2 



AT 



= o, 



R 

c/H 



T>(a far ) -D(a ncar ) 



"h 



(25) 



Using this we calculate the expected ISW signal for the 
50 high est significance voids in Table 4 of IGranett et al.l 
(2008a), employing the concordance ACDM cosmology to 
determine af ar and a near for each void from the void red- 
shift measurements. The ISW signal is found to be only 
—0.42 jjK on average if the dark matter underdensities are 
smaller than the observed underdensities in the LRG counts 
by the bias factor of 2.2 (taking as — 0.8). This is in con- 
trast to the detected mean signal of — 11.3pK which is over 
20 times bigger! We must therefore conclude that the void 
radii and/or underdenities have been significantly underes- 
timated. The void radii can at most be increased by a factor 
of 1.75 within the quoted uncertainties so the observed sig- 
nal of — 11.3 fiK can be matched only if the underdensities 
are increased by a factor of 5 (implying a bias factor of 0.2). 
The CMB temperature decrements of such model voids cal- 
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Figure 2. The variation with increasing void radius of the variance of the Hubble parameter, the density contrast, the density parameter 
and the peculiar velocity for the ACDM power-law and CHDM bump models, given the WMAP-5 and SDSS data (with la limits). 



culated using eq. (|25[) are shown in Fig[7] and are (by con- 
st ruction) similar to the a ctual measurements shown in Fig. 2 
of lGranett et al.l (|2008bh . While such an underbids for the 
observed LRGs may seem implausible, we emphasise that 
this is the o nly way in wh i ch the temperature decrements 
observed by iGranett et al.l ([2008a, b) can be accounted for 
as being due to the late ISW effect. 

FiglZl displays a histogram of the probabilities for find- 
ing such voids in the SDSS LRG survey volume (5 h~ 3 Gpc 3 
in the redshift range 0.4 < z < 0.75). The most improbable 
void is at z — 0.672 — in order to yield the observed av- 
erage CMB temperature decrement it must have a density 
contrast of -0.72 (quoted galaxy underdensity of -0.316 mul- 
tiplied by 5/2.2) and a radius of 230 h" 1 Mpc (radius derived 
from the quoted volume of 10 7 h~ s Mpc 3 and multiplied by 
1.75). The probability of such a void is 1.9 x 10~ 247 accord- 
ing to our calculations. Although linear theory may not be 
applicable for such a deep void, it is clear that its existence 
is in gross conflict with the standard theory of structure 
formation from gaussian primordial density perturbations. 

This conclusion is strengthened by the recent detection 
of very large peculiar velocities on large scales. As seen in 
FigOD the expected variance of the peculiar velocity as cal- 
culated by eq. (|14[) is about 200 km s _1 on a scale of 100 /i _1 
Mpc, whereas the measured value is at least 5 times higher, 



and t he discrepancy is even b igger on larger scales up to 300 
Mpc JKashlinskv et alj|2008h . 

It is also seen from Fig[5]that if a determination of the 
Hubble constant is required with say 1% accuracy, then mea- 
surements extending out to at least 150 h~ x Mpc must be 
made to overcome local fluctuations. A similar estimate was 
made bv lLi et al.l (|2008j) who noted that the observed vari- 
ance in measurements of h is in accord, thus consistent with 
the assumption of a gaussian density fi eld. However the voids 
observed in the SDSS LRG survey (JGranett et al.l l2008bl ) 
call this assumption into question. In particular whether 
there is a large local void is then an issue that must be 
addressed observationally and not dismissed on the grounds 
that it is inconsistent with gaussian perturbations. The Hub- 
ble flow is presently poorly measured in the redshift range 
0.1 < -s < 0.3 — just w here the effects of such a local void 
would be most apparent 11 Alexander et al.ll2007f ). Given that 
dark energy may well be an artifact of such a void, this issue 
needs urgent attention. 

The question of how such voids can have been gener- 
ated without conflicting with the CMB observations is be- 
yond the scope of the present work. Some suggestions have 
been made in the context of multi-field inflationary models 
dOcchionero et al.l I1997J ; iDiMarco fc Notaril I2006J ; lltzhaki I 
l200Sfl . 
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Figure 3. The probability distribution of the Hubble contrast (with \a limits), given the WMAP-5 and SDSS data, for the ACDM 
power-law and CHDM bump models, for spherical voids of radius R = (40, 70, 100, 150, 200, 300, 500, 800) xh,- 1 Mpc. 
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Figure 5. The probability of a fluctuation in the Hubble contrast greater than or equal to a given value lP H in a sphere of radius R (with 
Iff limits), given the WMAP-5 and SDSS data, for the ACDM power-law and CHDM bump models for R = (40, 70, 100, 150, 200, 300, 
500, 800) xr'Mpc. 
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Figure 6. The probability of a fluctuation in the density contrast less than or equal to a given value S° in a sphere of radius R (with 
la limits), given the WMAP-5 and SDSS data, for the ACDM power-law and CHDM bump models for R = (40, 70, 100, 150, 200, 300, 
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Figure 7. The left panel shows the ISW signals of the 50 voids detected bv lGranett et al.l I l2008af) . calculated using eq,l l25l l; in order 
to match the observed average ISW signal of — 11.3/iK it has been necessary to increase the void radii by a factor of 1.75 and the 
underdensities by a factor of 5. The right panel show the probability of such voids occurring in the SDSS LRG survey volume according 
to the concordance ACDM model. 
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